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Most working scientists hold fast to the concept of 'realism' - a viewpoint according to which 
an external reality exists independent of observation. But quantum physics has shattered some 
of our cornerstone beliefs. According to Bell's theorem, any theory that is based on the joint 
assumption of realism and locality (meaning that local events cannot be affected by actions in space- 
like separated regions) is at variance with certain quantum predictions. Experiments with entangled 
pairs of particles have amply confirmed these quantum predictions, thus rendering local realistic 
theories untenable. Maintaining realism as a fundamental concept would therefore necessitate the 
introduction of 'spooky' actions that defy locality. Here we show by both theory and experiment 
that a broad and rather reasonable class of such non-local realistic theories is incompatible with 
experimentally observable quantum correlations. In the experiment, we measure previously untested 
correlations between two entangled photons, and show that these correlations violate an inequality 
proposed by Leggett for non-local realistic theories. Our result suggests that giving up the concept of 
locality is not sufficient to be consistent with quantum experiments, unless certain intuitive features 
of realism are abandoned pQ. 

Physical realism suggests that the results of observations are a consequence of properties carried by physical 
systems. It remains surprising that this tenet is very little challenged, as its significance goes far beyond science. 
Quantum physics, however, questions this concept in a very deep way. To maintain a realistic description of 
nature, non-local hidden-variable theories are being discussed as a possible completion of quantum theory. They 
offer to explain intrinsic quantum phenomena - above all, quantum entanglement [2] - by nonlocal influences. 
Up to now, however, it has not been possible to test such theories in experiments. We present an inequality, 
similar in spirit to the seminal one given by Clauser, Home, Shimony and Holt [3^ on local hidden variables, that 
allows us to test an important class of non-local hidden-variable theories against quantum theory. The theories 
under test provide an explanation of all existing two qubit Bell- type experiments. Our derivation is based on a 
recent incompatibility theorem by Leggett [3], which we extend so as to make it applicable to real experimental 
situations and also to allow simultaneous tests of all local hidden-variable models. Finally, we perform an experiment 
that violates the new inequality and hence excludes for the first time a broad class of non-local hidden- variable theories. 

Quantum theory gives only probabilistic predictions for individual events. Can one go beyond this? Einstein's 
view [5, 6 was that quantum theory does not provide a complete description of physical reality: "While we have thus 
shown that the wavefunction does not provide a complete description of the physical reality, we left open the question 
of whether or not such a description exists. We believe, however, that such a theory is possible." [5]. It remained an 
open question whether the theory could be completed in Einstein's sense [7J. If so, more complete theories based on 
objective properties of physical systems should be possible. Such models are referred to as hidden-variable theories. 

Bell's theorem [B] proves that all hidden-variable theories based on the joint assumption of locality and realism 
are at variance with the predictions of quantum physics. Locality prohibits any influences between events in 
space-like separated regions, while realism claims that all measurement outcomes depend on pre-existing properties 
of objects that are independent of the measurement. The more refined versions of Bell's theorem by Clauser, Hornc, 
Shimony and Holt [3] and by Clauser and Home [HI EI] start from the assumptions of local realism and result in 
inequalities for a set of statistical correlations (expectation values), which must be satisfied by all local realistic 
hidden-variable theories. The inequalities are violated by quantum mechanical predictions. Greenberger, Home and 
Zcilinger [TTJ [12] showed that already perfect correlations of systems with at least three particles are inconsistent 
with these assumptions. So far, all experiments motivated by these theorems are in full agreement with quantum 
predictions E3 EH E3 EH EO EH] ■ For some time, loopholes existed that allowed the observed correlations to be 
explained within local realistic theories. In particular, an ideal Bell experiment has to be performed with detectors 
of sufficiently high efficiency (to close the 'detection loophole') and with experimental settings that are randomly 
chosen in space-like separated regions (to close the 'locality loophole'). Since the first successful Bell experiment 
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by Freedman and Clauser (T3], later implementations have continuously converged to closing both the locality 
loophole [15j [THJ [HE [20] on the one hand and the detection loophole [17l [21] on the other hand. Therefore it is 
reasonable to consider the violation of local realism a well established fact. 

The logical conclusion one can draw from the violation of local realism is that at least one of its assumptions 
fails. Specifically, cither locality or realism or both cannot provide a foundational basis for quantum theory. Each 
of the resulting possible positions has strong supporters and opponents in the scientific community. However, Bell's 
theorem is unbiased with respect to these views: on the basis of this theorem, one cannot, even in principle, favour 
one over the other. It is therefore important to ask whether incompatibility theorems similar to Bell's can be found 
in which at least one of these concepts is relaxed. Our work addresses a broad class of non-local hidden-variable 
theories that are based on a very plausible type of realism and that provide an explanation for all existing Bell-type 
experiments. Nevertheless we demonstrate, both in theory and experiment, their conflict with quantum predictions 
and observed measurement data. Following the recent approach of Leggett [4] , who introduced the class of non-local 
models and formulated an incompatibility theorem, we have analysed its assumptions and derived an inequality valid 
for such theories that can be experimentally tested. In addition, the experiments allow for a simultaneous test of all 
local hidden-variable models - that is, the measurement data can neither be explained by a local realistic model nor 
by the considered class of non-local models. 

The theories under investigation describe experiments on pairs of particles. It is sufficient for our purposes to 
discuss two-dimensional quantum systems. We will hence focus our description on the polarization degree of freedom 
of photons. The theories are based on the following assumptions: (1) all measurement outcomes are determined 
by pre-existing properties of particles independent of the measurement (realism); (2) physical states are statistical 
mixtures of subensembles with definite polarization, where (3) polarization is defined such that expectation values 
taken for each subensemble obey Malus' law (that is, the well-known cosine dependence of the intensity of a polarized 
beam after an ideal polarizer). 

These assumptions are in a way appealing, because they provide a natural explanation of quantum mechanically 
separable states (polarization states indeed obey Malus' law). In addition, they do not explicitly demand locality; that 
is, measurement outcomes may very well depend on parameters in space-like separated regions. As a consequence, such 
theories can explain important features of quantum mechanically entangled (non-separable) states of two particles 
(a specific model can be found in Appendix I): first, they do not allow information to be transmitted faster than 
the speed of light; second, they reproduce perfect correlations for all measurements in the same bases, which is a 
fundamental feature of the Bell singlet state; and third, they provide a model for all thus far performed experiments 
in which the Clauser, Home, Shimony and Holt (CHSH) inequality was violated. Nevertheless, we will show that all 
models based on assumptions (l)-(3) are at variance with other quantum predictions. 

A general framework of such models is the following: assumption (1) requires that an individual binary 
measurement outcome A for a polarization measurement along direction a (that is, whether a single photon is 
transmitted or absorbed by a polarizer set at a specific angle) is predetermined by some set of hidden-variables A, 
and a three-dimensional vector u, as well as by some set of other possibly non-local parameters 77 (for example, 
measurement settings in space-like separated regions) - that is, A — A(\,u,a,rj). According to assumption (3), 
particles with the same u but with different A build up subensembles of 'definite polarization' described by a 
probability distribution /05(A). The expectation value A{u), obtained by averaging over A, fulfils Malus' law, 
that is, A(u) = J dAps(A)^4(A, u, a, rf) = u ■ a. Finally, with assumption (2), the measured expectation value for 
a general physical state is given by averaging over the distribution F(u) of subensembles, that is, (A) — J duF{u)A{u). 

Let us consider a specific source, which emits pairs of photons with well-defined polarizations u and v to laboratories 
of Alice and Bob, respectively. The local polarization measurement outcomes A and B are fully determined by the 
polarization vector, by an additional set of hidden variables A specific to the source and by any set of parameters rj 
outside the source. For reasons of clarity, we choose an explicit non-local dependence of the outcomes on the settings a 
and b of the measurement devices. Note, however, that this is just an example of a possible non-local dependence, and 
that one can choose any other set out of 77. Each emitted pair is fully defined by the subensemble distribution Pu,v{^)- 
In agreement with assumption (3) we impose the following conditions on the predictions for local averages of such 
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measurements (all polarizations and measurement directions are represented as vectors on the Poincare sphere [22 ]): 

A(u) = J dXp s ^(X)A(a,b,X) = u- a (1) 
B(v) = / d\ Pa ^(X)B(b, a, A) =v-b (2) 



It is important to note that the validity of Malus' law imposes the non-signalling condition on the investigated 
non-local models, as the local expectation values do only depend on local parameters. The correlation function of 
measurement results for a source emitting well-polarized photons is defined as the average of the products of the 
individual measurement outcomes: 

AB(u,v) = J dXp s ^(X)A(a,b,X)B(b,a,X) (3) 

For a general source producing mixtures of polarized photons the observable correlations are averaged over a distri- 
bution of the polarizations F(u,v), and the general correlation function E is given by: 

E = (AB) = J dudvF(u, v)AB(u, v) (4) 

It is a very important trait of this model that there exist subensembles of definite polarizations (independent of 
measurements) and that the predictions for the subensembles agree with Malus' law. It is clear that other classes 
of non-local theories, possibly even fully compliant with all quantum mechanical predictions, might exist that do 
not have this property when reproducing entangled states. Such theories may, for example, include additional 
communication [23] or dimensions [23]. A specific case deserving comment is Bohm's theory [25]. There the non-local 
correlations are a consequence of the non-local quantum potential, which exerts suitable torque on the particles 
leading to experimental results compliant with quantum mechanics. In that theory, neither of the two particles in a 
maximally entangled state carries any angular momentum at all when emerging from the source 26J. In contrast, 
in the Leggett model, it is the total ensemble emitted by the source that carries no angular momentum, which is a 
consequence of averaging over the individual particles' well defined angular momenta (polarization). 

The theories described here are incompatible with quantum theory. The basic idea of the incompatibility theorem ^ 
uses the following identity, which holds for any numbers A = ±1 and B = ±1: 

- 1 + \A + B\ = AB = 1 - \A- B\ (5) 

One can apply this identity to the dichotomic measurement results A = A(a,b,X) = ±1 and B = B(b,a,X) = ±1. 
The identity holds even if the values of A and B mutually depend on each other. For example, the value of a specific 
outcome A can depend on the value of an actually obtained result B. In contrast, in the derivation of the CHSH 
inequality it is necessary to assume that A and B do not depend on each other. Therefore, any kind of non-local 
dependencies used in the present class of theories are allowed. Taking the average over the subensembles with definite 
polarizations we obtain: 

- 1 + J dXp S M(X)\A + B\ = J dXp^{X)AB = 1 - J dXp^{X)\A - B\ (6) 

Denoting these averages by bars, one arrives at the shorter expression: 

- 1 + \A + B\ =AB = 1 - \A-B\ (7) 
As the average of the modulus is greater than or equal to the modulus of the averages, one gets the set of inequalities: 

-1 + \A + B\ < AB<1-\A-B\ (8) 

By inserting Malus' law, equations and in equation and by using expression Q, one arrives at a set of 
inequalities for experimentally accessible correlation functions (for a detailed derivation see Appendix II). In particular, 
if we let Alice choose her observable from the set of two settings a\ and a^, and Bob from the set of three settings b%, 
&2 and &3 = 02, the following generalized Leggett-type inequality is obtained: 

Snlhv = |£nM + £ 23 (0)| + \E 2 2{v) + E 23 (0)\ < 4 - -| sin £| (9) 

7T 2 
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where Eki(<p) is a uniform average of all correlation functions, denned in the plane of 3k and bi, with the same relative 
angle ip; the subscript NLHV stands for 'non-local hidden-variables'. For the inequality to be applied, vectors a\ and 
bi necessarily have to lie in a plane orthogonal to the one defined by 0,2 and 62. This contrasts with the standard 
experimental configuration used to test the CHSH inequality, which is maximally violated for settings in one plane. 




FIG. 1: Testing non-local hidden-variable theories, a, Diagram of a standard two-photon experiment to test for hidden- 
variable theories. When pumping a nonlinear crystal (NL) with a strong pump field, photon pairs are created via SPDC and 
their polarization is detected with single-photon counters (PC). Local measurements at A and B are performed along directions 
a and b on the Poincare sphere, respectively. Depending on the measurement directions, the obtained correlations can be 
used to test Bell inequalities (b) or Leggett-type inequalities (c). b, Correlations in one plane. Shown are measurements 
along directions in the linear plane of the Poincare sphere (H (V) denotes horizontal (vertical) polarization). The original 
experiments by Wu and Shaknov [27] and Kocher and Commins [28], designed to test quantum predictions for correlated 
photon pairs, measured perfect correlations (solid lines). Measurements along the dashed line allow a Bell test, as was first 
performed by Freedman and Clauser [13]. c, Correlations in orthogonal planes. All current experimental tests to violate Bell's 
inequality (CHSH) are performed within the shaded plane. Out-of-plane measurements are required for a direct test of the 
class of non-local hidden-variable theories, as was first suggested by Leggett. 

The situation resembles in a way the status of the Einstein, Podolsky and Rosen (EPR) paradox before the advent 
of Bell's theorem and its first experimental tests. The experiments of Wu and Shaknov [37] and of Kocher and 
Commins [ 128] were designed to demonstrate the validity of a quantum description of photon-pair correlations. As this 
task only required the testing of correlations along the same polarization direction, their results could not provide 
experimental data for the newly derived Bell inequalities (Fig. [TJl, b). Curiously, as was shown by Clauser, Horne, 
Shimony and Holt, only a small modification of the measurement directions, such that non-perfect correlations 
of an entangled state are probed, was sufficient to test Bell's inequalities. The seminal experiment by Freedman 
and Clauser [13] was the first direct and successful test [29]. Today, all Bell tests - that is, tests of local realism - 
are performed by testing correlations of measurements along directions that lie in the same plane of the Poincare 
sphere. Similar to the previous case, violation of the Leggett-type inequality requires only small modifications to that 
arrangement: To test the inequality, correlations of measurements along two orthogonal planes have to be probed 
(Fig. [lb). Therefore the existing data of all Bell tests cannot be used to test the class of nonlocal theories considered 
here. 

Quantum theory violates inequality ([9]). Consider the quantum predictions for the polarization singlet state of two 
photons, \^>~)ab = ^75 ~ |^0a|-£0,b], where, for example, \H) a denotes a horizontally polarized photon 

propagating to Alice. The quantum correlation function for the measurements and bi performed on photons 
depends only on the relative angles between these vectors, and therefore E^i — — ■ b\ — — cos (p. Thus the left hand 
side of inequality ([9|, for quantum predictions, reads |2(cos<^ + 1)|. The maximal violation of inequality (|9| is for 
(fmax = 18.8°. For this difference angle, the bound given by inequality ^ equals 3.792 and the quantum value is 3.893. 

Although this excludes the non-local models, it might still be possible that the obtained correlations could be 
explained by a local realistic model. In order to avoid that, we have to exclude both local realistic and non-local realistic 
hidden- variable theories. Note however that such local realistic theories need not be constrained by assumptions (1)- 
(3). The violation of the CHSH inequality invalidates all local realistic models. If one takes 

Schsh = I + E12 — E 2 i + E22I < 2 (10) 
the quantum value of the left hand side for the settings used to maximally violate inequality S is 2.2156. 
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The correlation function determined in an actual experiment is typically reduced by a visibility factor V to 
E exp = — V cos (p owing to noise and imperfections. Thus to observe violations of inequality ^ (and inequality (10)) 
in the experiment, one must have a sufficiently high experimental visibility of the observed interference. For the 
optimal difference angle ip ma x — 18.8°, the minimum required visibility is given by the ratio of the bound (3.792) 
and the quantum value (3.893) of inequality ([9|, or ~ 97.4%. We note that in standard Bell-type experiments, a 
minimum visibility of only ~ 71% is sufficient to violate the CHSH inequality, inequality ( 10 1, at the optimal settings. 
For the settings used here, the critical visibility reads 2/2.2156 ~ 90.3%, which is much lower than 97.4%. 




Bob 



Pump laser 
>.=395 nm 



FIG. 2: Experimental set-up. A 2-mm-thick type-II /3-barium-borate (BBO) crystal is pumped with a pulsed frequency- 
doubled Tksapphire laser (180 fs) at A = 395 nm wavelength and ~ 150 mW optical c.w. power. The crystal is aligned to 
produce the polarization-entangled singlet state \^~)ab = b — \ V)a\H) b\- Spatial and temporal distinguishability 

of the produced photons (induced by birefringence in the BBO) are compensated by a combination of half-wave plates (A/2) 
and additional BBO crystals (BBO/2), while spectral distinguishability (due to the broad spectrum of the pulsed pump) is 
eliminated by narrow spectral filtering of 1 nm bandwidth in front of each detector. In addition, the reduced pump power 
diminishes higher-order SPDC emissions of multiple photon pairs. This allows us to achieve a two-photon visibility of about 
99%, which is well beyond the required threshold of 97.4%. The arrows in the Poincare spheres indicate the measurement 
settings of Alice's and Bob's polarizers for the maximal violation of inequality Note that setting 62 lies in the y-z plane and 
therefore a quarter-wave plate has to be introduced on Bob's side. The coloured planes indicate the measurement directions 
for various difference angles <p for both inequalities. 



In the experiment (see Fig. |2|, we generate pairs of polarization entangled photons via spontaneous parametric 
down-conversion (SPDC). The photon source is aligned to produce pairs in the polarization singlet state. We 
observed maximal coincidence count rates (per 10 s), in the H/V basis, of around 3,500 with single count rates 
of 95,000 (Alice) and 105,000 (Bob), 3,300 coincidences in the ±45° basis (75,000 singles at Alice and 90,000 at 
Bob), and 2,400 coincidences in the R/L basis (70,000 singles at Alice and 70,000 at Bob). The reduced count 
rates in the R/L basis are due to additional retarding elements in the beam path. The two-photon visibilities 
are approximately 99.0 ± 1.2% in the H/V basis, 99.2 ± 1.6% in the ±45° basis and 98.9 ± 1.7% in the R/L 
basis, which - to our knowledge - is the highest reported visibility for a pulsed SPDC scheme. So far, no exper- 
imental evidence against the rotational invariance of the singlet state exists. We therefore replace the rotation 
averaged correlation functions in inequality ^ with their values measured for one pair of settings (in the given plane). 

In terms of experimental count rates, the correlation function E(a, b) for a given pair of general measurement 
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settings is denned by 

E(a,b)= N ++ + N -- N +-- N -+ (11) 

where Nab denotes the number of coincident detection events between Alice's and Bob's measurements within the 
integration time. We ascribe the number +1, if Alice (Bob) detects a photon polarized along a (b), and —1 for the 

orthogonal direction a 1 - (ir~). For example, N^ denotes the number of coincidences in which Alice obtains a and 

Bob b- 1 . Note that E(ctk,bi) — Eki(<p), where ip is the difference angle between the vectors a and b on the Poincare 
sphere. 



To test inequality Q9J> , three correlation functions (E(ai,bi,), E(d 2 ,b 2 ), E{a 2 ,b 3 )) have to be extracted from the 
measured data. We choose observables a% and b\ as linear polarization measurements (in the x-z plane on the 
Poincare sphere; see Fig. 2| and 02 and b 2 as elliptical polarization measurements in the y-z plane. Two further 



correlation functions {E(a,2, bi) and £7(ai,&2)) are extracted to test the CHSH inequality, inequality (10 1 



The first set of correlations, in the x-z plane, is obtained by using linear polarizers set to ct\ and j3\ (relative to 
the 2-axis) at Alice's and Bob's location, respectively. In particular, a\ — ±45°, while j3\ is chosen to lie between 
45° and 160° (green arrows in Fig. [2|. The second set of correlations (necessary for CHSH) is obtained in the same 
plane for a 2 = 0°/90° and (3\ between 45° and 160°. The set of correlations for measurements in the y-z plane is 
obtained by introducing a quarter-wave plate with the fast axis aligned along the (horizontal) 0°-direction at Bob's 
site, which effectively rotates the polarization state by 90° around the z-axis on the Poincare sphere (red arrows 
in Fig. [2]). The polarizer angles are then set to 02 = 0°/90° and [3 2 is scanned between 0° and 115°. With the 
same (3 2 and a\ = ±45°, the expectation values specific only for the CHSH case are measured. The remaining 
measurement for inequality ([9| is the check of perfect correlations, for which we choose 012 = P3 = 0°, that is, 
the intersection of the two orthogonal planes. Figure [3] shows the experimental violation of inequalities (|9| and 
( |10| for various difference angles. Maximum violation of inequality ^ is achieved, for example, for the settings 
{«!, a 2 , /3x, 02, 13 3 } = {45°, 0°, 55°, 10°, 0°}. 




FIG. 3: Experimental violation of the inequalities for non-local hidden-variable theories (NLHV) and for local 
realistic theories (CHSH). a, Dashed line indicates the bound of inequality (JoJ for the investigated class of nonlocal hidden- 
variable theories (see text). The solid line is the quantum theoretical prediction reduced by the experimental visibility. The 
shown experimental data were taken for various difference angles ip (on the Poincare sphere) of local measurement settings. 
The bound is clearly violated for 4° < ip < 36°. Maximum violation is observed for ipmax ~ 20°. b, At the same time, no 
local realistic theory can model the correlations for the investigated settings as the same set of data also violates the CHSH 
inequality ( ] 1 Op . The bound (dashed line) is overcome for all values p around p ma x, and hence excludes any local realistic 
explanation of the observed correlations in a. Again, the solid line is the quantum prediction for the observed experimental 
visibility. Error bars indicate s.d. 

We finally obtain the following expectation values for a difference angle tp = 20° (the errors are calculated assuming 
that the counts follow a poissonian distribution): E(ai,h) = -0.9298 ± 0.0105, E(a 2 ,b 2 ) = -0.942 ± 0.0112, 
E(a 2l b 3 ) = -0.9902 ± 0.0118. This results in S NLH v = 3.8521 ± 0.0227, which violates inequality {9} by 
3.2 standard deviations (see Fig. [3]). At the same time, we can extract the additional correlation func- 
tions E(a 2 ,b 1 ) = 0.3436 ± 0.0088, E(a!,b 2 ) = 0.0374 ± 0.0091 required for the CHSH inequality. We obtain 
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Schsh = 2.178 ± 0.0199, which is a violation by ~ 9 standard deviations. The stronger violation of inequality (10 1 
is due to the relaxed visibility requirements on the probed entangled state. 

We have experimentally excluded a class of important non-local hidden-variable theories. In an attempt to 
model quantum correlations of entangled states, the theories under consideration assume realism, a source emitting 
classical mixtures of polarized particles (for which Malus' law is valid) and arbitrary non-local dependencies via the 
measurement devices. Besides their natural assumptions, the main appealing feature of these theories is that they 
allow us both to model perfect correlations of entangled states and to explain all existing Bell-type experiments. We 
believe that the experimental exclusion of this particular class indicates that any non-local extension of quantum 
theory has to be highly counterintuitive. For example, the concept of ensembles of particles carrying definite 
polarization could fail. Furthermore, one could consider the breakdown of other assumptions that are implicit in our 
reasoning leading to the inequality. These include Aristotelian logic, counterf actual definiteness, absence of actions 
into the past or a world that is not completely deterministic |30j . We believe that our results lend strong support to 
the view that any future extension of quantum theory that is in agreement with experiments must abandon certain 
features of realistic descriptions. 
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APPENDIX I: AN EXPLICIT NON-LOCAL HIDDEN- VARIABLE MODEL 

We construct an explicit non-local model compliant with the introduced assumptions (l)-(3). It perfectly simulates 
all quantum mechanical predictions for measurements in a plane of the Poincare sphere. We model the correlation 
function of the singlet state, E^ 1 — —a-b, for which all local averages (A) and (B) vanish. In particular, the violation 
of any CHSH-type inequality can be explained within the model and, in addition, all perfect correlations state can be 
recovered. 

Let us start with a source that emits photons with well-defined polarization. Polarization u is sent to Alice and v 
to Bob. Alice sets her measuring device to a and Bob to b. The hidden-variable A € [0, 1] is carried by both particles 
and predetermines the individual measurement result as follows: 

A = A(a, u, A) = {+; £ A A ||° A ^j; with A A = 1(1 +H-0), (12) 

where A is the outcome of Alice. This means, whenever A < Xa the result of the measurement A is +1, and for 
A > Xa the result is —1. Note that the measurement settings only enter in Xa and are hence independent of the 
hidden- variable A of the source. The outcome of Bob is given by 

B = B(a,b,u,v, A) = { ^ , n with x\, x% € [0, ll arbitrary but xo — x\ = — (1 + v- &T13) 

v ' ' ' ' ' I -1 for A <E [0, xi) U (x 2 , 1], J 2 

All non-local dependencies are put on the side of Bob. His measuring device has the information about the setting 
of Alice, a, and her polarization u. The requirement of the non-local models discussed here is that the local averages 
performed on the subensemble of definite (but arbitrary) polarizations u and v obey Malus' Law, i.e. A^ = u ■ a for 
Alice, and B$ — v ■ b for Bob. Indeed, a straight-forward calculation shows that this requirement is fulfilled for both 
Alice and Bob: 

A u = I dX- I dX = 2X A -l=u-a, (14) 

JO J\ A 
r%2 r^i /•! 

B$= I dX- dX- dX = 2{x 2 -x 1 )-l = v-b. (15) 

Thus the current construction fulfills one of our aims - we recover Malus' law. In order to get the correct formula for 
correlated counts one can fix the value of x\ and x 2 in the following way: 

X\ = -[l + u-a — v-b + a-b], 

x 2 = - [3 + u ■ a + v ■ b + a ■ b]. (16) 

With this definitions whenever x\ < Xa < x 2 the expectation value for measurements on the subensembles reproduces 
quantum correlations. Simply: 

rxi i-Xa rX2 pi ^ 

A S B$ = - dX+ dX- dX+ dX = 2(X A - a?i - x 2 ) + 1 = -a ■ b. (17) 

^0 J Xi J \a J X2 

Therefore, in the next step, one must find the conditions for which both x\ and x 2 take values from [0, 1] and X\ and 
x 2 take values from [0, 1]. 



To this end, using definitions ( 16 1 one finds that the first condition is equivalent to a set of four inequalities: 



— 1 + u • 6 < a ■ b + u ■ a < 3 + v ■ b, 

-3-v-b< a-b + u-a < 1 -v-b. (18) 

Note that the upper bound 3 + v ■ b cannot be exceeded by the middle term, as well as the lower bound — 3 — v ■ b. 
Thus, this set of four inequalities is equivalent to a single one: 

\a-b + u-a\<l-v-b. (19) 
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Similarly, the second condition can be reexpressed as: 

\a-b-u-a\<l + v-b. (20) 



Finally, the validity condition for the model is a conjunction of (19 1 and (20 1 



\a-b±u-a\<lTv-b. (21) 

If this relation is not satisfied the model does not recover quantum correlations. Either it becomes inconsistent since 
xi or x 2 leave their range or the necessary relation x\ < Xa < x 2 is not satisfied, or both. This is the origin of the 
incompatibility with general quantum predictions. Nevertheless the model can explain all perfect correlations and the 
violation of the CHSH inequality. 

Imagine a source producing pairs with the following property: whenever polarization u is sent to Alice polarization 
v = — u is sent to Bob. Both parties locally observe random polarizations. For Alice, the local average over different 
polarizations yields 

(A) = ^As+ ^Als = \u-a~ \u-a = 0, (22) 

as it should be for the singlet state. The same result holds for Bob. In this way, we have reproduced the randomness 
of local measurement outcomes, typical for measurements on entangled states. 

With the same source, one can explain perfect correlations for measurements along the same basis, i.e. b — ±o. To 
see how the model works take v = —u and b = a, and find that X\ = | [1 + u ■ a] = Xa and x 2 = 1. As it should be, 
Bob's outcomes are always opposite to Alice's: 

B = B(a,a,u,-u,X) = ( ^ for ^^^y (23) 

If in the same subensemble we take b = —a we obtain x\ = and x 2 — Xa, which results in B — A, again in 
full agreement with quantum mechanics. Note that for these measurement settings, i.e. in order to obtain perfect 



correlations, condition (21 1 imposes no additional restrictions since it is always satisfied. For example, if u is sent to 
Alice and b = —a one obtains \ — l±u-a\<l^fu-a, which always holds. The same argument applies to the other 
subensemble and other measurement possibilities b = ±a. 

Finally, the full predictions of quantum theory are recovered if Alice and Bob restrict their mea sure ments to lie 



in planes orthogonal to vectors u and v, respectively, i.e. u ■ a = v ■ b = 0. In this case, condition (211 is satisfied, 



as \a ■ b\ < 1. In general, if condition (21) is satisfied, i.e. for a consistent set of parameters, our model reproduces 
quantum correlations since they are already reproduced in every subensemble and hence averaging over different 
polarizations is not necessary: 



(AB) = A S B_ S = A_ a B a = -a-b. (24) 
Therefore every experimental violation of the CHSH inequality can be explained by the presented non-local model. 



APPENDIX II: DERIVATION OF THE INEQUALITY 

After Leggett 0], one can take a source which distributes pairs of well-polarized photons. Different pairs can have 
different polarizations. The size of a subensemble in which photons have polarizations u and v is described by the 
weight function F{u,v). All polarizations and measurement directions are represented as vectors on the Poincare 
sphere. In every such subensemble individual measurement outcomes are determined by hidden variables A. The 
hidden variables are distributed according to the distribution ps,v(X). 

For any dichotomic measurement results, A — ±1 and B — ±1, the following identity holds: 

-1 + \A + B\=AB = 1-\A-B\. (25) 

If the signs of A and B are the same | A + B\ = 2 and \A-B\=0, and if A = -B then \A + B\ = and \A—B\ = 2. 
Any kind of non-local dependencies is allowed, i.e. A = A(a, b, u, v, X, ....) and B = B(a, b, u, v, A, ....). Taking the 
average over the subensemble with definite polarizations gives: 



1 + J dX Pa ^(X)\A + B\ = J dXpa, ff {X)AB = 1 - J dX P ^(X)\A - B\, 



(26) 
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which in an abbreviated notation, where the averages are denoted by bars, is 



-1 + \A + B\ = AB = 1-\A-B\. (27) 

Since the average of the modulus is greater or equal to the modulus of the averages one gets the set of inequalities 

-1 + \A + B\ < AB<1-\A-B\. (28) 

From now on only the upper bound will be considered, however all steps apply to the lower bound as well. 
With the assumption that photons with well defined polarization obey Malus' law: 

A = u-a, B = v-b, (29) 



the upper bound of Eq. (28 1 becomes: 

AB<l-\H-S k -v-bi\, (30) 

where 3k and bi are unit vectors associated with the fcth measurement setting of Alice and the Zth of Bob, respectively. 
Taking the average over arbitrary polarizations one obtains 

Ef-i<l- J dudvF(u,v)\u ■ a k - v -bi\, (31) 

where E k i is the correlation function which can be experimentally measured when Alice chooses to measure 3k and 
Bob chooses 6;. Let us denote by Uki and Vki the length of projections of vectors u and v onto the plane spanned by 
3k and bi . Since one can decompose vectors u and v into a vector orthogonal to the plane of the settings and a vector 
within the plane the scalar products read: 

u-3 k = u k icos((j) ak - <f> u ), (32) 
v-bi = v k i cos(<j) bl - cj> v ), (33) 

where all the <p angles are relative to some axis within the plane of the settings; Angles 4> u and <f> v describe the position 
of the projections of vectors u and v, respectively, whereas angles <p ak and 05 ; describe the position of the setting 
vectors. With this notation the inequality transforms to: 



E H < I - J dudvF(u, v)\uu cos(0 Qfe - cf> u ) - v k i cos(4> bl - 4> v )\. (34) 

The magnitudes of the projections can always be decomposed into the sum and the difference of two real numbers: 

um = rii + 7i 2 , vu = n-i - n 2 . (35) 
We insert this decomposition into the last inequality, and hence the terms multiplied by ri\ and ni read: 

cos(0 Qfc - (j) u ) - cos(0 bi - <j) v ) 

_ ■ + <t>bi - {<l>u + <j>v) . ~M + <K - <Pv ,oa\ 

— z sin ^ sm 2 ' v / 

cos(0 Qfc - <j) u ) + cos(0 6i - 4>v) 

= 2 cos + *h - fa + M cos ^ ~ ^ -(^-M 

2 2 y ' 

respectively. We make the following substitution for the measurement angles: 

4fei = ^ ' Vkl = ^ ak " ( S8 > 

and parameterize the position of the projections within their plane by: 

Ipuv = S , Xuv = Pu-<Pv (39) 
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Using these new angles one arrives at: 

Eki(^khfki) < 1 - 2 / dudvF(u, v)\n 2 cos ^ - Xuv cos(£ fc ; - ip uv ) - m sin LPkl Xuv sin(£ fei - ip uv )\, (40) 



where in the correlation function Eki(£ki, Pki) we explicitly state the angles it is dependent on. The expression within 
the modulus is a linear combination of two harmonic functions of — ip uv , and therefore is a harmonic function 



itself. Its amplitude reads \ n\ cos 2 ( yfci 2 X "" ) + n 2 sin 2 Xtt " ), and the phase is some fixed real number a 



E kl ^ k i,ip k i) <l-2 / dudvF(ti,v)^n 2 2 co S 2 ( ^ - Xuv ) + n?sin 2 ( ^ ^ )[ cos(^ - V™ + <*)!• (41) 

In the next step we average both sides of this inequality over the measurement angle £ki = ^ ak ^ bl . This means an 
integration over £ki € [0, 2tt) and a multiplication by The integral of the £m dependent part of the right-hand side 



of (41 1 reads: 



o 

By denoting the average of the correlation function over the angle £ki as 

(•27T 

2tt 



1 f 27r 2 

— / d&i | cos(£ fe i - <ip uv + a) | = -. (42) 



Eki(<Pkl) = ^- / d£,kiEki(£kh(fki), (43) 



o 



one can write (41 ) as 



< 1 - ^ / dAW-FW 3) ^1 cos2(^^) + n? sin 2 (^^). (44) 

This inequality is valid for any choice of observables in the plane defined by dk and bi. One can introduce two new 
observable vectors in this plane and write the inequality for the averaged correlation function Ek'i'iv'k'i') °f these new 
observables. The sum of these two inequalities is 

E k i(ipki) + E kn ,(ip' k , v ) < 2 - J J dudvF(u,v) 

x (y rii ccs3^ (45) 
One can use the triangle inequality 

\W+y\\ < 11*11 + 11*11, (46) 

^{x 1 +y 1 y 2 + {x 2 + y 2 y < yjx$ + xl + yjy^ + yl (47) 
for the two-dimensional vectors x — {x 1,2:2) and y = (y 1,2/2), with components defined by: 

I 'Phi ~ Xuv I I fk'l' ~ Xuv I /Ao\ 

2:i = |n 2 cos 1, y x = |n 2 cos K 1, (48) 

I • fhl ~ Xuv I I . '■Ph'V ~ Xuv I r An \ 

x 2 = \n l sm 1, y 2 = \ni sin 1. (49) 

This implies that the integrand is bounded from below by: 



n 2 C0S 2 m - Xuv + n2 s . n 2 m - Xuv + ^2 cQg2 gfc^ _ Xgg + ^ ^2 ^ _ Xuv 



>\n$(\cos o \ + \cos^— 1 +nj | sin 1 + | sin 1 . (50) 
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One can further estimate this bound by using the following relations 



|cos( )| + |cos( )| > | sin 1 and (51) 

I • / ^fikl Xuv \ I | . ff'k'l' ~ Xuv,, . | . Phi ^Pk'l' I /co \ 

|sm( )| + |sm(-^— )| > | sin — (52) 



This estimate follows if one uses the formula for the sine of the difference angle to the right-hand side argument 

^Z*Vii = m-x.. _ fU'-x.. Nam el y , 



Isin ^'-^' l = | sin ^ ~ Xuv cos ^ - X ™ - cos ^' - Xuv sin - Xuv | (53) 
2 2 2 2 2 

^ I • <y5fci — Xuv n f'k'l'~Xuv, . I Vfei — Xuv n ■ f'k'l'~Xuv, /■_,-. 

< I sin 1| cos 1 + I cos ||sin-^ 1. (54) 



After these estimates, the lower bound of E k i + Ek'i 1 (following form the left-hand side inequality in (28)) is equal 
to minus the upper bound, and thus one can can apply the upper bound to the modulus of the left hand side of 
( 45 1 . This is because the only formal difference between expressions in the estimates seeking the lower bound of the 
averaged Eq. (28 1 compared to those seeking the upper bound boils down to the interchange between ri\ and TI2. 
After applying (51 I and (52), this makes no difference anymore. One can shortly write: 

\Ekl(Vkl) + Ek'l'(<p' k n')\ < 2- ^|sin ^ fc '~^ fc ''' | J dudvF(u,v)^/n 2 + n\. (55) 
Going back to the magnitudes: 

\E u {y u ) + E wv Wm>)\ < 2-^\sm ipkl ~ 2 P ' k ' l ' \ J dudvF(u,v)^/u 2 kl + v 2 kl . (56) 

This inequality is valid for any choice of the plane of observables. The bound involves only the projections of vectors 
u and v onto the plane of the settings. The integrations in the bound can be thought of as a mean value of expression 
V u ki + v ki averaged over the distribution of the vectors. For the plane orthogonal to the initial one the inequality is 

\E^ pq ) + E^ q ,(^ plq ,)\ < 2- 2 ^\sm Vpq ~ 2 p/p ' q ' \ f dudvF{u,v)^u 2 pq + v 2 pq . 



(57) 

where u pq and v pq denote the projections of vectors u and v, respectively, onto the plane spanned by the settings 
a p and b q (which is by construction orthogonal to the plane spanned by a k and b{). We add the inequalities for 
orthogonal observation planes, (56 1 and (57 1, and choose (p' k , v = tp' p i q i = and ifki — <p pq = (p. This gives 

\E U {V) + E kn >m + \E pq W) + E^ q ,(0)\ < 4 - ^| sin || J dudvF(u,v) (*fifa+tf> l + yjuj q + v 2 ^ . (58) 



We apply the triangle inequality (46 1 to the expression within the bracket. This time vectors x and y have the 
following components: 

x*=(u k i,u pq ), y= (v k i , Vpq). (59) 

The integrand is bounded by: 

\l U ll + V il + \] U lq + '"lq>\J { u kl + U p q) 2 + (v kl + V pq ) 2 (60) 

Let us consider the term involving vector u only. Since the lengths are positive 

(uki+u pq ) 2 >u 2 kl + u 2 pq . (61) 

Recall that Uki and u pq are projections onto orthogonal planes. One can introduce normal vectors to these planes, 
Tiki and n pqi respectively, and write 

(rtki ■ u) 2 + u 2 kl = 1, and (n pq ■ u) 2 + u 2 = 1. (62) 
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Note that the scalar products are two components of vector u in the Cartesian frame build out of vectors flki, n pq , 
and the one which is orthogonal to these two. Since vector u is normalized one has: 

(n fei • u) 2 + (n pq ■ u) 2 < 1, (63) 

which implies for the sum of equations ( |62[ ) : 

u 2 l+ u 2 pq >l. (64) 
The same applies to vector v and one can conclude that 



u li + v li + V< + < > V^- (65) 



Since the weight function F(u,v) is normalized, the final inequality reads: 



\E kl &)+E kn ,(0)\ + \E^)+E± ql (Q)\<4-^\ S m^\. (66) 



APPENDIX III: CORRIGENDUM TO THE PREVIOUS VERSION ON QUANT-PH (0704.2529V1) 

The experimental values of correlation functions were measured at the difference angle ip = 20° (on the Poincare 
sphere) rather than 18.8° as stated in the previous version. The reported violation of 3.2 standard deviations refers 
to ip = 20°. 

In Appendix II the boundaries of all integrals over variables x an d "0 should be the following: \ varies from — 27r 
to 2ir; tp varies from |x|/2 to 2ir — |x|/2. Also, in equation (46) and following F(8 U ,9 V ) should be replaced by 

lXdxJ^ xU2 d^F{e u ,9 v ^,x). 

The final inequality of the paper and all its conclusions do not change. The errors arose in the attempt for a 
concise presentation but were never used in the actual computation. In the present manuscript these errors are 
corrected. We replaced Appendix II with a shorter derivation of the final inequality. 

We would like to thank Stephen Parrott for pointing us to the errors in the original paper. 
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